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Abstract 

The dynamics of the one-dimensional spin- 1/2 quantum XXZ model with random fields is inves- 
tigated by the recurrence relations method. When the fields satisfy the bimodal distribution, the 
system shows a crossover between a collective-mode behavior and a central-peak one with increas- 
ing field while the anisotropy parameter A is small; a disordered behavior replaces the crossover 
as A increases. For the cases of Gaussian and double-Gaussian distributions, when the standard 
deviation is small, the results are similar to that of the bimodal distribution; while the standard 
deviation is large enough, the system only shows a disordered behavior regardless of A. 
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1. INTRODUCTION 



Quantum spin systems are of considerable interest for the reasons that they provide a 
ground for studying quantum many-particle phenomena, and offer the possibility to compare 
the theoretical and experimental results, etc. The dynamical properties of the quantum spin 
systems have received much attention and a variety of achievements have been got over the 
years. Among the systems studied, the one-dimensional (1-D) spin-1/2 XXZ chain is one of 
the nontrivial models and_has been used to describe several quasi-l-D compounds such as 



CsCoBr 3 



lrfl, CsCoCl 3 J El |, Cs 2 CoCl 4 



, |6( and TiCoCl3 [7]. Concerning the dynamics 



of this model, the spin correlation function has been studied by using approximation and 
numerical methods 8Ml7|| . Strong numerical evidence was found for a change of the bulk- 
spin autocorrelation function at T = 00 from Gaussian decay to exponential decay for the 
anisotropy parameter A increasing from 0, and from exponential decay to power-law d ecay 



as A approaches 1, then replaced by a more rapid decay upon further increase of A 



id- 



The dynamics of the equivalent-neighbor XXZ model was studied in much detail at T = 
and T = 00 using different calculational techniques, and was found the same long-time 
asymptotic behavior for the correlation function [16l . Il7j ] . 

More efforts have been concentrated on the dynamics of random quantum spin systems 



in the past decades which can be used 
electric crystals 



18 



21] and spin glasses 



or describing more real materials such as ferro- 



22J. Florencio and Barreto have studied the ran- 



dom transverse Ising model and obtained that the system undergoes a crossover between a 
collective-mode behavior and a central-peak one while the exchange couplings or external 
fields satisfy the bimodal distribution [23] . Later, the dynamics of the four-body transverse 



Ising model has been investigated for the cases of bond and field randomness [24], following 
which, a series of spin systems such as the XY model with zero magnetic, the transverse Ising 
model with Gaussian distribution and two-dimensional transverse Ising model etc. have also 
been studied 25M29] . A recent study shows that the next-nearest-neighbor interaction has a 
strong influence on the dynamics of the Ising system 30] . For the random XXZ model, the 
spin correlations have been inve stig ated by using exact diagonalization 31], the real space 



renormalization group method 



32 



34] and a finite-chain study [35l . 136] etc. Infinite tem- 



perature spin-spin correlation function has been found to display exponential localization 
in space indicating insulating behavior for large enough random fields 32]. The transverse 



2 



correlation function at T = has been found to exhibit a power-law decay to exponential 
decay depending on the exchange disorder jsg]. 

In this paper, we investigate the effects of the random fields on the time evolution of 
the quantum XXZ model in the high-temperature limit. We find that the system with the 
random fields that satisfy the bimodal distribution undergoes a crossover between a central- 
peak behavior and a collective-mode one with increasing field when the anisotropy parameter 
A is small (e.g., A = 0.01), but the collective-mode behavior vanishes as A approaches 0.4, 
then when A increases to 1.0, the central-peak behavior vanishes and the system just shows 
a disordered behavior. 

The arrangement of this paper is organized as follows. In Sec. |2] we give a simple 
introduction to the 1-D quantum XXZ model and the recurrence relations method. In Sec. 
|3]we discuss the results, and Sec. @] contains a summary. 



2. MODEL AND METHOD 

The Hamiltonian of the 1-D quantum spin- 1/2 XXZ model with external fields can be 
written as 

i i 

where af (a = x, y, z) are Pauli spin operators, J and A are the exchange coupling and the 
anisotropy parameter, respectively. Bi denote the external fields, which may be regarded as 
random variables. Clearly, this Hamiltonian can describe two special cases: the Ising model 
for which A = oo and the isotropy XY model where A = 0. 

The spin autocorrelation function plays an important part in the study of the dynamics 
of quantum spin systems. It is defined as 

C {t) = <rr; (/) rr; ((»)> (2) 

where (• • ■ ) denotes an ensemble average followed by an average over the random variable. 
The corresponding spectral density which is the Fourier transform of C (t) can be expressed 
as 

$ (u) = / e wt C (t) dt, (3) 
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and for the mathematical simplicity, $ (u) is able to be obtained as 



$ (u) = lim 

£->0 



POO 

Re / dtC(t)e 
Jo 



zt 



(4) 



where z = e + iu, e > 0. 

The recurrence relations method has 



23 



jeen proved to be very powerful in the calculation of 



24 



27H30I. l37f . Next, we will give a brief introduction 



dynamic correlation function 13 
to this method. 

Considering a Hermitian operator a? (t) as a dynamical variable and then expanding it 
with an orthogonal set in a Hilbert space, 



"5(*) = 5>" 



(5) 



u=0 



where a u (t) are the time-dependent coefficients. There exists the following set of recurrence 
relation for the basis vectors f v , 



f u+1 = iLf v + A„f v -i, v>0, 

\fvi fv) 



A,. 



V > 1, 



{fv-l, fv-l) 

where L = [H, ] is the quantum Liouvillian operator, (f u , f u ) - 
A = 1. The coefficients a u (t) in Eq. (jSJ) satisfy the relation: 

A,y +1 a u+1 (t) = -a u [t) + a^-i (t) , v > 0, 

where a v (t) = — — — , a_i (t) = 0. 



(6) 
(7) 



and 



(8) 



The spin autocorrelation function C (t) can be expressed as the form of moment expansion 



with 



,2k 



fc=0 v ; 



1 



,r =-Tra*[H, [H, ■ ■ ■ [H, ^] • • •]] : 



(9) 



(10) 



where fi 2k is the 2fcth moment of C (t). Supposing that the first Q moments have been 
calculated by Eq. ( ITU]) , we can obtain C (t) by constructing the pade approximate. Because 
of the mathematical complexity, a finite number of moments can be got, the large times of 
C (t) are divergent, thus we just can discuss the short times of C (t). 



By taking the inner product for (t) and crj (0), from Eq. (j3j), we can find that do (t) is 
the spin autocorrelation function 



a (t) = (o?(t)o?(0)> = C(t). 

Applying the Laplace transform a v (z) = J °° e~ zt a u (t) dt [z = e + iu,e > 0) to Eq.® , a (z) 
in the continued-fraction representation can be obtained as 

a (z) = ^ . (11) 

z + k~ 

z + 



z + • • • 

Because we can only get a finite number of recurrants, it is necessary to terminate it with 



some schemes. Here, we use the Gaussian terminator |13|, |38| Which can best serve our 
problem. With the help of the Pade approximate and the Gaussian terminator, we can 
obtain the spin autocorrelation function C (t) and the corresponding spectral density $ (cj) 
of the system, respectively. 



3. RESULTS AND DISCUSSIONS 



In order to investigate the spin autocorrelation function C (t) for crj, we choose the zeroth 
basis vector /o = a? . With the recurrence relation Eq. ([6]), the remaining basis vectors can 
be obtained: 

fx = Bjrf + ./AnV , - Jo) ,o ; - Ja* j+1 a* + JAa»o* +1 , 

f 2 = (A, - B) - 2 J 2 - 2J 2 A 2 ) a] + U 2 Ao]_ x + 2J 2 AaJ +1 + .l 1 Ao ] a) y. 

t2 a „x JU i t2 „x „y jy o t2 „x jy „y , t2 „x „y „y 

-J Ao i _ 2 a j _ x a, J + J a 1+1 a j _ 1 a j - 2J a^^a^ + J a^a^ 

-J 2 Aa* +2 ay j+l + J 2 Aa* +l ay j+2 - 2B j jAa*a*_ l + {!',, . + B,) Jo*.^ 

+ (B^ + 5,) Ja^a) + J 2 AaJ_ l( x^_ 2 - J 2 ^.^^ + J 2 Aa-| +1 cr|cr|_ 1 

-25, JA^ +1 - 2./ 2 AV;Vr: ; rr;.. + J 2 A^_ lf r|^ +1 - •/ V. .ajrrj . , 

+J 2 Aa* +1 cx|cx J 2 +2 , 

etc. The first three norms of the basis vectors are obtained as follows: 

(/o,/o) = l, 
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(A, A) = B] + 2 J 2 + 2 J 2 A 2 , 



(A, A) = A 2 - 2AiBj + - 4Ai J 2 + B 2 _ V J 2 + 2B j - 1 B j J 2 + 65 2 J 2 + 2B j+1 BjJ 2 
+B 2 +1 J 2 + 12 J 1 _ 4A X J 2 A 2 + 12B 2 J 2 A 2 + 24 + 8 7W. 

Then the continued fraction coefficients can be got from Eq. (j7]). 

Next, numerical results of the spin autocorrelation functions C (t) and the spectral den- 
sities $ (w) are given when the external fields satisfy three types of distributions: bimodal 
distribution, Gaussian distribution and double-Gaussian distribution. With special values 
of the anisotropy parameter A, the effects of the random external fields on the dynamics of 
the given system are investigated as follows. 

3.1. Bimodal distribution 

We first consider the case that the external fields satisfy the bimodal distribution. 

P (Bi) = p5 (Bi - B x ) + (l-p)S (B { - B 2 ) . (12) 

For simplicity and without loss of generality, we choose the exchange coupling J = 1.0, 
which sets the energy scale, and the external fields B\ = 1.8, B 2 = 0.2. For A = 0.01, 0.1, 
0.4 and 1.0, the results of the spin autocorrelation function C (t) and the spectral density 
$ (u) are shown in Fig. 1 and Fig. 2, respectively. The continued-fraction coefficients are 
presented in the insets. 

When A is small (e.g., A = 0.01, 0.1), the spin autocorrelation function 
[see in Figs. l(al), (a2)], changes from a monotonically decreasing behavior to a damped 
oscillatory one as p increases. When p = 0, the exchange coupling energy is higher than the 
external field energy. The interaction among the spins is stronger than that between the 
spin and the external field. The dynamics is dominated by the exchange coupling energy 
and shows a central-peak behavior. When p = 0.25, the spin autocorrelation function has a 
slight fluctuation and the fluctuation becomes acute with the increase of the external fields. 
When p = 1, the value of the external field is larger than that of the exchange coupling. 
The system behaves as the precession of independent spins about the field and the exchange 
coupling causes a damping. Hence, the system presents a collective-mode behavior. Figs. 
2(bl) and (b2) show that the peak of $ (u) moves from oj = to 2 as p increases, which also 

6 



reaches the conclusion that the system undergoes a crossover from the central-peak behavior 
to the collective-mode one with increasing field when A is small. 

Figure l(a3) shows that when A = 0.4, the dynamics of the system changes from a 
central-peak regime to a disordered behavior which is intervenient between a central-peak 
one and a collective-mode one as p increases. The spin autocorrelation function for p = 
decays monotonically to and the spectral density is now peaked at u = 0. So, the system 
is at the central-peak regime where the dynamics is mostly dominated by the exchange 
coupling. By comparing the curve for p = in Figs. l(al) or (a2) to the one in Fig. I(a3), 
we find that C (t) for A = 0.4 decays faster than that for A = 0.01 or A = 0.1. As p 
increases to 1.0, the system shows not a collective- mode behavior but a disordered one. The 
spectral density displayed in Fig. 2(b3) tends to have an expansion at high frequency. 

When A = 1.0 [see Figs. I(a4) and 2(b4)], the system presents a disordered behavior as 
the concentration of B\ increases, i.e., in this case, the dynamics of the system can not be 
characterized by either behavior singly. Specially the case of p = is very similar to the 



most-disordered case mentioned in Ref. 23]. By comparing the results when A = 0.01 to 
that of the 1-D XY model, we find that they are very similar, so the effect of the anisotropy 
parameter can be basically ignored, in other words, the dynamics of the system is governed 
by the competition between spin-spin interactions and the external fields. Comparing Figs. 
l(al), (a2), (a3) and (a4), it can be found that as the anisotropy parameter A increases, the 
crossover from the central-peak behavior to the collective-mode one vanishes. When A = 
1.0, the system becomes the 1-D quantum Heisenberg system. The anisotropy parameter 
together with the external field and the exchange coupling decide the dynamic behaviors 
of the system. The competition between the spin-spin interactions and the external fields 
becomes very fierce which drives the system to be disordered. 

3.2. Gaussian distribution 

In the case, the external fields satisfy the Gaussian distribution, 

P (B. t ) = -=L- exp [- {B t - Bf /2a%] , (13) 



where B is the mean value of the external fields, Ob is the standard deviation. Here, we 
take the exchange coupling J equal to 1.0 and B equal to 0.0, 0.5, 1.0, 1.5, 2.0. When the 
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anisotropy parameter A = 0.1, 0.4 and 1.0, the results of the spin autocorrelation function 
and the spectral density are displayed in Fig. 3 for <jb = 0.3 and Fig. 4 for ob = 3.0. 

Figure 3(al) shows that when the standard deviation is small (ctb = 0.3) and A = 0.1, 
the system shows two types of dynamics as B increases: the central-peak behavior and the 
collective-mode behavior. In this case, the effect of A can be ignored basically, the dynamics 
of the system changes according to the concentration of B. As A increases from 0.1 to 1.0, 
a disordered behavior replaces the central-peak behavior or the collective-mode behavior. 
Figs. 3(a2), (b2) and (c2) also show that the effects of the external fields are to urge the 
system to show a crossover when A = 0.1 and as A increases to 1.0, the system displays a 
disordered behavior. 

When the standard deviation is large enough (cr B = 3.0) [see Fig. 4], the system just 
shows a disordered behavior which is something in between the central-peak behavior and 
the collective-mode one regardless of the anisotropy parameter. From Fig. 3 and Fig. 4, it 
is not difficult to see that both the crossover and disordered behavior are replaced by one 
disordered behavior with the increase of ob- This is because the value of the external field 
is large and the value range is wide when ob = 3.0. The large external fields drive the spin 
orientation of the system to be disordered. 

3.3. Double- Gaussian distribution 

Double-Gaussian distribution is a common form of bimodal distribution and Gaussian 
distribution, which can be used to describe both discrete distribution and continuous one, 

P (B { ) = P^L- exp [- {Bi - B,) 2 /2a 2 B ] + (1 - p) - ? =L- exp [- (B { - B 2 f /2a*] , 

(14) 

where < p < 1 represents the concentration of B\ that satisfies the Gaussian distribution. 
The external fields satisfy Eq. ( TT4"|) . in which the mean values B\ = 1.8 and B 2 = 0.2, 
and the exchange coupling is constant (J = 1.0). C (t) and $ (u) for A = 0.1, 0.4 and 1.0 
are calculated and the results are shown in Fig. 5 and Fig. 6, respectively. The figures 
indicate that when ob = 0.3, the system for A = 0.1 shows a crossover between a central- 
peak behavior and a collective-mode one, and a disordered behavior as A increases to 1.0. 
However, the system only shows a disordered behavior when <Jb = 3.0, no matter what A 
takes. 
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From above discussion, we can see that the dynamical behavior of the system is affected 
by the competition between the spin-spin interactions and the external fields, not by the 
different disordered distribution. Also, it is easy to find that the dynamics of the system is 



28 



when A is small (e.g., A = O.Of). When 



similar to that of the 1-D quantum XY model 
A = 0, the XXZ model becomes the isotropy XY model. We find that the above results are 
the same as those in Ref. 



28 



when we take A = 0. 



4. SUMMARY 



We have studied the dynamics of the 1-D spin-1/2 quantum XXZ model in the random 
external fields at the high-temperature limit by means of the recurrence relations method. 
We find that the dynamics of the system with three types of random distributions are 
affected by the competition among the external field, the anisotropy parameter and the 
exchange coupling, but the anisotropy parameter can be basically ignored when it is small 
(e.g., A = 0.01). For the case of bimodal disorder, when the anisotropy parameter A is small 
(e.g., A = 0.01), the dynamics of the system undergoes a crossover between a collective-mode 
behavior and a central-peak one with increasing field; as A increases to 0.4, the dynamics of 
the system changes from a central-peak regime to a disordered behavior which is intervenient 
between a central-peak one and a collective-mode one with the increase of the fields; then as 
A approaches 1, the system shows a disordered behavior. In the cases of Gaussian disorder 
and double-Gaussian disorder, when the standard deviation of the random field <7b is small, 
a disordered behavior replaces the crossover as A increases. When ob becomes large enough, 
the system shows only a most disordered behavior regardless of the anisotropy parameter. 
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Figure Captions 



Fig. 1 The spin autocorrelation functions where the external fields take the values B\ = 
1.8 with probability p and B 2 = 0.2 with probability (1 — p) . (al), (a2), (a3), (a4) correspond 
to the cases of A = 0.01, 0.1, 0.4 and 1.0. The continued-fraction coefficients are presented 
in the insets. 

Fig. 2 The corresponding spectral densities for the same parameters as in Fig. 1. (bl), 
(b2), (b3), (b4) correspond to the cases of A = 0.01, 0.1, 0.4 and 1.0. 

Fig. 3 The spin autocorrelation functions and the spectral densities for the case that the 
external fields satisfy the Gaussian distribution. The standard deviation <tb takes 0.3. (al), 
(a2) correspond to the case that A equals to 0.1, (bl), (b2) correspond to the case that A 
equals to 0.4, and (cl), (c2) correspond to the case that A equals to 1.0. 

Fig. 4 The spin autocorrelation functions and the spectral densities for the case that the 
external fields satisfy the Gaussian distribution. The standard deviation Ob takes 3.0. (al), 
(a2) correspond to the case that A equals to 0.1, (bl), (b2) correspond to the case that A 
equals to 0.4, and (cl), (c2) correspond to the case that A equals to 1.0. 

Fig. 5 The spin autocorrelation functions for the case that the external fields satisfy 
the double-Gaussian distribution in which the mean values B\ = 1.8 and B 2 = 0.2 with 
probabilities p and (1 — p). (al), (a2) and (a3) correspond to the case that the standard 
deviation a B takes 0.3 and A takes the value 0.1, 0.4, 1.0, respectively, (bl), (b2) and (b3) 
correspond to the case that the standard deviation Ob takes 3.0 and A takes the value 0.1, 
0.4, 1.0, respectively. 

Fig. 6 The spectral densities for the same parameters as in Fig. 5. (al), (a2) and (a3) 
correspond to the case that the standard deviation Ob takes 0.3 and A takes the value 0.1, 
0.4, 1.0, respectively, (bl), (b2) and (b3) correspond to the case that the standard deviation 
<j b takes 3.0 and A takes the value 0.1, 0.4, 1.0, respectively. 
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